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ABSTRACT 

We present absorption-line index gradients for a sample of SO galaxies in the Fornax 
Cluster. The sample has been selected to span a wide range in galaxy mass, and the 
deep VLT-F0RS2 spectroscopy allows us to explore the stellar populations all the 
way to the outer disk-dominated regions of these galaxies. We find that globally, in 
both bulges and disks, star formation ceased earliest in the most massive systems, 
as a further manifestation of downsizing. However, within many galaxies, we find an 
age gradient which indicates that star formation ended first in the outermost regions. 
Metallicity gradients, when detected, are always negative such that the galaxy centres 
are more metal-rich. This finding fits with a picture in which star formation continued 
in the central regions, with enriched material, after it had stopped in the outskirts. 
Age and metallicity gradients are correlated, suggesting that large differences in star 
formation history between the inner and outer parts of SO galaxies yield large differ- 
ences in their chemical enrichment. In agreement with previous results, we conclude 
that the radial variations in the stellar populations of SO galaxies are compatible with 
the hypothesis that these galaxies are the descendants of spiral galaxies whose star 
formation has ceased. With the addition of radial gradient information, we are able to 
show that this shutdown of star formation occurred from the outside inward, with the 
later star formation in the central regions offering a plausible mechanism for enhancing 
the bulge light in these systems, as the transformation to more bulge-dominated SO 
galaxies requires. 

Key words: galaxies: elliptical and lenticular - galaxies: evolution - galaxies: forma- 
tion - galaxies: stellar content - galaxies: clusters: individual: NGC 1316:NGC 1380 



1 INTRODUCTION 

In order to better understand the formation and evolution 
of SO galaxies, we have undertaken a detailed study of these 
systems in the Fornax Cluster, combining deep optical long- 
slit spectroscopy obtained with the FORS instrument on the 
VLT with archival optical and near-infrared imaging. The 
sample and major-axis kinematics of these galaxies are pre- 
sented in Bedregal et al. (2006a, hereafter Paper I). In Bedre- 
gal et al. (2006b, hereafter Paper II), we used those data to 
study the TuUy-Fisher relation (TuUy & Fisher 1977) for 
these systems as part of a larger compilation of 60 local 
SOs, and concluded that their disk kinematics were consis- 
tent with what would be expected if these systems were 
simply faded SO galaxies. In Bedregal et al. (2008, hereafter 
Paper III) we studied the central stellar population prop- 
erties of these galaxies, and found evidence that the bulge- 
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dominated regions of these galaxies had star-formation prop- 
erties that varied systematically with the mass of their host 
galaxies. In this paper, we seek to tie together these bulge 
and disk properties by studying the gradients in common 
absorption-line indices. By obtaining such gradients out as 
far as the region where the disk dominates the light, and in- 
terpreting them in terms of simple stellar population models, 
we seek an integrated understanding of the formation and 
evolution of both of the major components of SO galaxies. 

It is already well known that elliptical galaxies and SO 
systems show variations in the properties of their stellar pop- 
ulations with radius. Early studies using broad-band colours 
established that the nuclear regions are usually redder than 
the outskirts, suggesting radial variations in their stellar 
population properties (Kormendy & Djorgovski 1989, and 
references therein). The first study of gradients in spectral 
absorption features was performed by McClure (1969): in 
a sample of 7 galaxies, he found that an index measuring 
the central CNA4216 band strength was stronger than at 
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larger radii, which he interpreted as the result of metallicity 
variations between the difTerent regions. Subsequent spectro- 
scopic studies, many based on the Lick/IDS line strength in- 
dices (Worthey et al. 1994), have explored absorption feature 
gradients using a broad range of line indices (e.g. Spinrad et 
al. 1972; Welch & Forrester 1972; Joly & Andrillat 1973; 
Cohen 1979; Couture & Hardy 1988; Peletier 1989; Gorgas 
et al. 1990; Bender & Surma 1992; Davies et al. 1993; Car- 
oUo et al. 1993; Gonzalez 1993; Fisher et al. 1995; Gorgas et 
al. 1997; Mehlert et al. 2003; Sanchez-Blazquez et al. 2006, 
2007). Overall, these studies indicate the presence of strong 
gradients in CN features (at 3883 and 4216 A), weaker gradi- 
ents in the Mg-triplet (at ~ 5200 A), g-band, Na D lines, Ca 
H&K features and some Fel sensitive indices, and no gra- 
dients in the hydrogen Balmer series (e.g. H/?, H7), MgH, 
TiO, Cal and the infrared Ca triplet. Many of these studies 
suggest that the existence of gradients in metallic features 
is produced by radial variations of metallicity, decreasing 
outward (e.g. Cohen 1979; Davies et al. 1993; Kobayashi & 
Arimoto 1999; Mehlert et al. 2003). Other works, however, 
suggest that there is also a gradient in the stellar popula- 
tion age, the central regions being younger than the outer 
parts (e.g. Gorgas et al. 1990; Munn 1992; Gonzalez 1993; 
Gonzalez & Gorgas 1996). 

These absorption-line index gradients provide clues to 
the mechanisms responsible of galaxy formation. Chemo- 
dynamical models of scenarios such as monolithic collapse, 
merger-driven evolution and other secular processes pre- 
dict different gradients. Although quantitative modeling of 
these effects is still in its infancy, clear qualitative indicators 
are becoming established. For example, monolithic collapse 
models (e.g. Eggen, Lynden-Bell & Sandage 1962; Larson 
1974a; Carlberg 1984; Arimoto & Yoshii 1987) generically 
predict the formation of strong gradients for both age- and 
metallicity-sensitive indices, while, in comparison, a forma- 
tion history driven by mergers will weaken any such gra- 
dients significantly (e.g. White 1980; Bekki & Shioya 1999; 
Barnes & Hernquist 1991; Mihos & Hernquist 1994). 

Clues to the more detailed star formation history of 
a galaxy can be obtained by comparing individual chemical 
elements. The different timescales associated with the super- 
nova producing a-elements (such as Mg) and Fe mean that 
two galaxies with the same mean stellar age but different 
histories will display different gradients in these elements. 
For example, a galaxy in which the stars were produced by 
a short starburst in the distant past and one where the star 
formation was over a more extended period may display sim- 
ilar constant Fe index values with radius, but the former will 
show a decreasing gradient in Mg line strengths while the 
latter will display a flatter Mg index profile (e.g. Thomas 
1999). 

Thus, absorption- line index gradients offer a valuable 
tool in the study of galaxy formation and evolution. Al- 
though we are not yet in a position to evaluate such gradi- 
ents in an unequivocal quantitative manner, their basic in- 
terpretation is fairly clear. Even without detailed modeling, 
their differential use is robust in that systematic differences 
between galaxies of different types are generally unambigu- 
ous. In addition, significant progress is being made with the 
quantitative analysis of such data, through the development 
of more sophisticated stellar population models. Hence the 
philosophy of this work is initially to make use just of the 



directly-observable line-strength gradients of the Fornax SOs 
to seek qualitative insight into any systematic variations in 
their evolutionary histories, and then to go on to use the 
new stellar population models of Vazdekis et al. (2010) to 
interpret such results in a more quantitative manner. 

The remainder of the paper is as laid out follows. In Sec- 
tion 2, we present the line index gradients derived from the 
deep long-slit spectroscopy of 9 Fornax Cluster SO galax- 
ies, with the associated detailed plots and values for indi- 
vidual galaxies provided in Appendix A. Section 3 describes 
the use of the new stellar population models to obtain ages 
and chemical abundances (and their gradients) from the 
observed line indices, with the detailed data given in Ap- 
pendix B. In Section 4, we discuss the implications of the 
results on the formation history of SO galaxies, while Sec- 
tion 5 summarizes these results. 



2 LINE STRENGTH INDEX GRADIENTS 

In this section, the index gradients measured for ten common 
absorption lines are presented. These are the same indices 
derived for the centres of these galaxies in Paper III, and the 
details of how they and their associated uncertainties were 
measured are presented in that paper. A description of how 
the spectra were binned in the spatial direction is also pre- 
sented in Paper III, Sec. 2. The binning was designed to yield 
a minimum signal-to-noise ratio in each spatial bin so that 
reliable spectral indices could be measured. See Paper HI for 
details. Note that for consistency all indices are calculated 
on a magnitude scale. 

In Appendix A, Figs. lAll to IA5I show line index mea- 
surements along the photometric major axes of the sample 
galaxies. The radial profiles have been folded around their 
galax;ies' centres to give a sense of any systematic variations. 
The radii have been scaled to each galaxy's half-light radius 
(hereafter, Rhalf) as calculated from infrared photometry 
(in the Ks-band; see Paper II). The straight lines show linear 
fits of these logarithmic parameters, with those points lying 
at radii small enough to be affected by seeing excluded from 
the fit. The resulting values of the index gradients and their 
uncertainties are presented in Table lXT] 

The line index gradients sample different fractions of 
the bulge- and disk-dominated regions of these galaxies. In 
Paper I we calculated a limiting radius (Rlim, dot-dashed 
blue lines in Figs. lAll to IA5|) from which the light pro- 
file is dominated by the disk component. Therefore, with 
Rlim we can assess the contribution of each component to 
the measured gradients. The relative bulge and disk con- 
tribution to the gradients is rather heterogeneous among 
our galaxy sample. For three galaxies the number of radial 
bins dominated by disk light is less than 30% (NGC 1380, 
1381 and 1375). For two other systems the disk contri- 
bution is between 70 and 90% (NGC 1380A and IC1963), 
while other two galaxies show gradients totally dominated 
by the disk light (ESO358-G006 and ESO 359-G002). Fi- 
nally, NGC 1316 and ESO 358-G059 have no information in 
this regard as Rlim cannot be estimated in those systems 
where the disk never dominates the luminosity profiles. 

In order to summarize this mass of information and 
search for systematic variations between Fornax Cluster SOs, 
Fig.[T] shows over-plotted the radial variation in all the line 
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Figure 1. Measurements of 10 line-strength indices as a function of radius for the entire galaxy sample. The key to the symbol representing 
each galaxy is shown at the top of the figure. The color code split the galaxy sample by luminosity into a bright (black), intermediate 
(blue) and faint (red) sample. For clarity, the error bars on each point, shown in the individual profiles of Figs. lAl1 - |A5l are not shown. 
Instead, in the top-right corner of each panel we present the median ±1 cr uncertainty of the measured indices. 
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Figure 2. Example of NGC 1316 radial data with Vazdekis et al. (2010) SSP model grids. H/3 is used as the age-sensitive index, while 
[MgFe]' (loft), (Fo) (center) and Mgb (right) are used as [M/H]-sensitive indicators. Color code on galaxy data represent the radial 
distance of each data point. The black asterisks represent the central measurements (within Rhalf/8). For clarity, only two model grids 
are shown, those with the higher (red) and lower (black) galaxy velocity dispersion. 



indices for all the sample galaxies. The galaxies have been 
divided and colour-coded by luminosity: 

• Bright SOs (-19.0 < Mb < -20.6): NGC1316, 1380 
and 1381, shown as black symbols. 

• Intermediate SOs (-18.1 Mb ^ -18.5): NGC 1380A, 
1375 and IC 1963, shown as blue symbols. 

• Famt SOs (-17.3 ^ Mb sS -17.5): ESO358-G059, 
ESQ 358-G006 and ESQ 359-G002, shown as red symbols. 

When interpreting the results we need to take into consid- 
eration that NGC 1316 (in the bright SO category) shows 
clear evidence of having suffered a recent merger (Schweizer 
1980; Goudfrooij et al. 2001a, b; Paper I), so may not be 
representative of the class. As is clear from Fig. [1] there are 
strong variations in behaviour both between different line 
indices and between different galaxies, but some systematic 
trends are immediately apparent. We therefore seek to ex- 
plore those systematic effects more fully. 

One initial point to note is that most of the stronger 
gradients measured are driven by the inner parts of bulge- 
dominated regions (e.g. Mg indices in NGC 1380, 1381, H/? 
in NGC 1375) . Generally, the radial profiles become fairly 
flat even before reaching IRhalf, immediately suggesting 
a scenario in which bulge and disk components may have 
undergone rather different mass assembly and star forma- 
tion histories. However, the radial behaviour of the indices 
present strong variations from galaxy to galaxy, suggesting 
a degree of variation in any such history. In addition, the dif- 
ferent line indices behave in rather different ways: all panels 
in Fig. [1] have vertical widths of 0.3 mag to allow a fair vi- 
sual comparison between radial index profiles, and so it is 
apparent that the H/3 and Mg indices display a wide variety 
of radial gradients, while the Fe tracers have much flatter 
profiles, but with a range of absolute values. 

Perhaps the most interesting phenomenon is that the 
H/3 index, while varying widely from galaxy to galaxy at 
small radii, seems to converge to a common low value at large 
radii. The lower-mass systems have systematically steeper 
gradients in this index, and correspondingly higher central 



values, while the higher-mass galaxies have similar low val- 
ues of H/3 at all radii. The simplest interpretation of this 
effect is that the lower-mass galaxies have undergone their 
final episode of star formation more recently, and that this 
star formation occurred primarily at small radii. It is pos- 
sible that the final star formation event in the higher-mass 
galaxies was similarly centrally concentrated, but occurred 
sufficiently long ago that any gradient in H/3 has now dis- 
appeared. We also need to be a little careful because vari- 
ations in the H/3 line index arise from both metallicity and 
age effects. We therefore now turn to interpreting the phys- 
ical significance of the H/3 profile and those of the other line 
indices. 



3 SINGLE STELLAR POPULATION 

PARAMETERS AND THEIR GRADIENTS 

Given the complex dependency of the observed absorption 
line indices on stellar population properties such as star- 
formation history and the chemical abundances of different 
elements, we rely on somewhat simplified models to trans- 
late the observations into parameters that have at least a 
crude physical interpretation. In particular, by fitting these 
indices to models of a single stellar population (SSP) with 
all stars having a common age and chemical abundances, we 
can derive a luminosity- weighted measure of the age, metal- 
licity and a-abundance of the observed population. 

In this section we will present the steps required to 
estimate these parameters using the new SSP models of 
Vazdekis et al. (2010) based on the MILES stellar library 
(Sanchez-Blazquez et al. 2006; Cenarro et al. 2007). This li- 
brary, with its ^1000 stellar spectra, provides very good 
coverage of the stellar atmosphere parameter space (grav- 
ity, temperature and metallicity) at a much higher spectral 
resolution than the original Lick/IDS stellar library: with an 
instrumental dispersion of « 56kms~^, it is well matched to 
the data analyzed here, allowing us, in many cases, to use the 
full observed spectral resolution when comparing data and 
models. Indeed, the velocity dispersion of the galaxy data 
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implies that in many situations we will be degrading the 
SSP models to the resolution of the galaxy spectra, and not 
the other way round, thus minimizing the loss of informa- 
tion. This optimal use of the data is particularly important 
in the fainter low-mass galaxies and in the disk-dominated 
outer regions of all the systems, as these regions share the 
observationally-challenging properties of low surface bright- 
ness and low intrinsic velocity dispersion. One further bene- 
fit of the MILES data is that they are flux calibrated, avoiding 
the complex and somewhat uncertain corrections involved in 
placing the absorption-line indices in the Lick/IDS library 
system. Note that the Mgi and Mg2 indices presented in Pa- 
per III were corrected to the Lick/IDS system, but the ones 
used here are not. 

The Vazdekis et al. (2010) models make full use of the 
MILES library and incorporate the latest developments in 
stellar evolution and SSP model techniques. In addition, 
they offer a publicly-available easily-usable web-based code 
(http://miles.iac.es/), which means that the analysis here 
can readily be reproduced and compared to the results from 
other data sets, thus making it the ideal choice for this anal- 
ysis. 

3.1 Stellar Population Parameters 

We estimate SSP ages and metallicities ([M/H]) using 
Vazdekis et al. (2010) model predictions and assuming a 
Kroupa-Revised initial mass function (IMF; Kroupa et 
al. 2003). The model spectra, corresponding to different com- 
binations of age and [M/H], were matched to the spectral 
resolution for each galaxy and radial position by convolving 
them with a Gaussian of dispersion 



gal 



(1) 



where Ugai is the velocity dispersion of the galaxy at the 
corresponding radius (see Paper I), and aub is the resolu- 
tion from the stellar library. Since aub has a small, but 
significant, wavelength dependence, the model spectra were 
degraded separately for each of the absorption lines used 
as the main diagnostics (H/3, Mg&, Fe5270 and Fe5335) 
using the value appropriate for the central wavelength of 
each feature. In the cases where cTgai < crub (NGC 1380A, 
IC 1963, ESQ 358-G059, ESQ 358-G006 and the outskirts of 
NGC 1381, NGC 1375 and ESQ 359-G002), the galaxy spec- 
tra were degraded to reach the library's resolution by con- 
volving them with a Gaussian of dispersion 



(2) 



where ainst is the F0RS2 instrumental dispersion of 
~30 kms-^ 

Once the model and galaxy spectra were matched in 
resolution, their line strength indices were measured. To 
estimate ages and [M/H], model grids for each individual 
galaxy data point were created as illustrated by the exam- 
ple shown in Fig. (2] In this figure we show the radial data of 
NGC 1316 with Vazdekis et al. (2010) SSP model grids for 
3 different metallicity tracers. The color code of the galaxy 
data represents the corresponding galactocentric distance. 
Note that we only show the grids for the models with the 
largest and smallest velocity dispersions to illustrate the am- 
plitude of this effect; for the actual fitting, each data point 




Figure 3. Gradients in Age, [M/H] and [a/Fe] tracer versus 
log(iJ;j -f^ax)' ^ proxy for dynamical mass, (see Paper III). Con- 
tinuous lines correspond to best linear fits. Dashed lines represent 
flat gradients. Symbols as in Fig.[l] 



was matched to its own grid of models. In the final modeling, 
the parameters were calculated using H/3 as the main age- 
sensitive index. Given the range of a-element abundances 
hinted at by the analysis of Section^ we adopt the com- 
bined [MgFe]' index as the primary measure of [M/H], due to 
its negligible sensitivity to a-element abundance (Gonzalez 
1993; Thomas, Maraston & Bender 2003). The best esti- 
mate SSP values of age and M/H] were then obtained from 
this grid using a bivariate polynomial interpolation between 
the calculated grid values. Finally, any variation with ra- 
dius in the relative Q-element abundances is traced using the 
Mg6/(Fe) ratio (see Paper III and Fig. 4 of Thomas, Maras- 
ton & Bender 2003). The resulting profiles for age, [M/H] 
and a-enhancement with radius for all of the sample galaxies 
are presented in Appendix B, Figs. lBll to lBSI 
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3.2 Uncertainties 

When interpreting any correlations apparent in tliese re- 
sults, it is important to realize that the nature of the model- 
ing process means that the estimated parameters are not in- 
dependent. We must therefore take some care in estimating 
the associated uncertainties. The errors in individual data 
points were calculated by Monte Carlo simulation, generat- 
ing pairs of values around the measured H/9 and [MgFe]' in- 
dices, with a dispersion matched to the observational errors 
on these quantities, placing them on the model grid, and 
then extracting the age and [M/H] parameters. The 68th 
percentile of the resulting age and [M/H] distributions then 
provides a realistic quasi-elliptical 1 a error contour around 
the derived Age- [M/H] value. The associated (often asym- 
metric) errorbars can then be derived by projecting the error 
contours onto the age and [M/H] axes. 

Similarly, the uncertainties in the Mg6/{Fe) ratio (our 
Q-abundance tracer) were estimated by Monte Carlo simu- 
lations, producing symmetric error bars in this case. 

3.3 Gradients 

Figs. lBlllBSI in Appendix B show the radial variation of the 
derived ages, [M/H] and a-enhancement indicators for all 
the galaxies, together with their uncertainties. Linear fits to 
those gradients are plotted as solid lines. The slopes of these 
fits (hereafter called "gradients") are given in Table iBl] 

The uncertainties in age and [M/H] gradients were cal- 
culated using the Monte Carlo simulations described above 
for every H/3-[MgFe]' pair, translated into simulated age- 
[M/H] pairs. From these values we created a large number 
of realizations of the radial distributions of these two pa- 
rameters for every galaxy, each of which was fitted with a 
linear function to give an estimate of the uncertainty in each 
gradient (and its correlation with other parameters). The 
quoted errors in these two parameters are asymmetric as a 
consequence of using individual age and [M/H] predictions 
having asymmetric uncertainties (see Sec. l3.2p . Errors in the 
a-enhancement gradients, also calculated via Monte Carlo 
simulations, are symmetric since no models are involved in 
their determination. 



4 RESULTS 

There are almost no previous studies of the stellar popula- 
tion gradients in SOs where these galaxies are treated as a 
separate class from the ellipticals (hereafter, Es). The work 
of Fisher, Franx & lUingworth (1996) is one of the few ex- 
ceptions. For a sample of 9 bright edge-on SOs they found 
that their bulge and disk metallicity gradients differ substan- 
tially, being steep and decreasing outwards in the bulge re- 
gion while becoming flat in the stellar disk. They also found 
that the central Mg and Fe indices correlate with central 
velocity dispersion (hereafter, ao) as in luminous Es. Hence, 
they conclude that SO bulges are probably the result of a dis- 
sipative collapse process (such as gas-rich mergers) at very 
early times. 

For a combined sample of 35 E and SO galaxies in 
the Coma cluster, Mehlert et al. (2003) calculate gradients 
reaching ~ IRhalf. They found that Es have on average 



slightly higher velocity dispersions, lower H/3 and higher 
metallic line-strengths than SOs. The latter form two fami- 
lies, one comparable to the Es and a second one with much 
younger stellar populations (~ 2 Gyr) , and weaker metal- 
lic lines. This result was later confirmed by Kuntschner 
et al. (2006) within the SAURON collaboration (Bacon et 
al. 2001). These authors studied a sample of 48 local Es and 
SOs and identified the young, low-cro SOs as the main drivers 
of the observed trends between gradients and central line 
indices. 

In the light of these results we now analyse the gradi- 
ents found in our Fornax SO sample. The first point apparent 
from Table [BTI is that virtually all the age gradients are pos- 
itive or fiat, indicating that the star formation ended first 
in the disk-dominated outer regions of the galaxies. This re- 
sult confirms the qualitative indication that we found from 
the H/3 gradients shown in Fig.[T] Indications of this re- 
sult were also found in Paper II, and we argued there that 
centrally-concentrated star formation at later times, enhanc- 
ing the bulge-to-disk ratio, is probably an integral stage in 
the morphological transformation of a spiral into a SO galaxy 
(Christlein & Zabludoff 2004; Shioya, Bekki & Couch 2004). 
Other observational evidence in favour of this scenario has 
been presented in the past (Moss & Whittle 2000; Mehlert 
et al. 2003; Poggianti et al. 2001; Bamford, Milvang- Jensen, 
& Aragon-Salamanca 2007). What is still unclear is whether 
this late central star formation represents an enhancement 
over the "quiescent spiral" star formation rate (perhaps fu- 
eled by interaction-induced inward gas fiows) or is just the 
last remaining centrally-concentrated star-formation gasp. 
The second hypothesis is probably favoured by some recent 
results (e.g.. Wolf et al. 2009) but the evidence is still incon- 
clusive. 

Figure[T] also suggested that the the cessation of star 
formation happened at earlier times in the brightest galax- 
ies, considerably reducing the H/3 and age gradients. If that 
trend is a reflection of an underlying galaxy mass depen- 
dency we would expect an anti-correlation between age gra- 
dient and galaxy mass. To test for this effect, Fig.[3]shows the 
various gradients plotted against Ra x Umaxj which provides 
a measure of each galaxy's mass (see Paper HI). There does 
seem to be a trend between age gradient and mass in the 
sense described (a Spearman rank test gives a probability 
of correlation between 95 and 97.5%), although NGC 1316, 
whose central regions have probably been rejuvenated by a 
recent merger, seems to deviate somewhat. Possible corre- 
lations between the gradients in [M/H] or a-enhancement 
and galaxy mass are either weaker or not detected (Spear- 
man rank tests give <90% and 90-95% probabilities, re- 
spectively). Previous studies (Mehlert et al. 2003; Rawle et 
al. 2010) found no correlations between these parameters 
in early-type galaxies (although see Spolaor et al. 2009 and 
Kuntschner et al. 2010 at low masses). In any case, some 
care is necessary in not over-interpreting these trends (or 
lack thereof) given the limitations imposed by the sample 
size and uncertainties. 

Close inspection of Figs. lBll to IB5I also reveals that the 
generally old age that we had ascribed to the disks in Sec- 
tion [2] is not the whole story, in that the age in the outer 
disk-dominated regions is > 10 Gyr for the brightest and 
most massive galaxies, while it is only ~ 4-5 Gyr for the less 
massive ones. This would seem to be a new manifestation 
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of "downsizing": star formation stopped first in tfie most 
massive galaxies not only in their centres [as found before in 
other early-type galaxy samples; see Smith at al. (2009) and 
references therein] but also in their oldest outermost regions. 

In Paper III we reported a tentative correlation between 
the a-element overabundance and age in the central regions 
of these galaxies. This correlation is confirmed here with the 
improved models. A Spearman rank test gives a probability 
of correlation between 95 and 97.5% (see Fig. |4] top-central 
panel). The correlation can be explained by considering a 
simple model where all galaxies started forming stars at ap- 
proximately the same time, formed stars for different time 
intervals At, and then ceased their star formation t Gyr ago. 
Since the SSP ages are correlated with t, while the a-element 
overabundance is anti-correlated with At, it follows that the 
ages must be correlated with the a-element overabundance, 
as observed (see Paper III and Aragon-Salamanca 2008 for 
details). Fig. |4] also confirms Paper III finding that there is 
no correlation between central ages and central [M/H], nor 
between central [M/H] and central a-enhancements. 

When detected, the [M/H] gradients are always neg- 
ative, implying that the stars in the inner regions of the 
galax;ies are on average more metal-enriched than those in 
the outskirts. This finding fits in with the growing picture 
that star formation continued, with enriched material, in the 
central regions after it had stopped in the outskirts. 

Gradients in the a-element enhancement tracer are usu- 
ally fiat within the uncertainties except for the two old and 
bright SOs (NGC 1380, 1381). These galaxies show clear neg- 
ative gradients (Table [Blj) . In the bottom-central panel of 
Fig.|4]we see a trend between Grad(Age) and Grad(a). This 
trend is consistent with the similar correlation found for the 
central values (a Spearman rank test gives a 90-95% proba- 
bility of correlation). 

Finally, the top-left panel of Fig.|4] clearly shows that 
gradients in age and [M/H] are correlated for the six bright- 
est galaxies (black and blue points). Nothing can be said for 
the faintest three, since no gradients were detected in them 
within the large errors (see Table IB1|I . A Spearman rank 
tests gives a 99.5% probability of correlation for all 9 SOs. 

It could be argued that, at least in part, this trend is 
driven by correlated errors as the 1 a error ellipses in Fig.|4] 
are clearly aligned with this trend. We have tested whether 
this correlation can be totally explained by the observed 
uncertainties using lO'* Monte Carlo simulations which con- 
sider the observed data points and their error distributions. 
Using a Kolmogorov-Smirnov test with a — 0.05 confidence 
level we find that 93% of the simulations reject the null hy- 
pothesis that our original sample and each simulated one 
were drawn from the same distribution. This indicates that, 
although there is a sizable correlation in the errors, there is 
a real correlation between Grad[M/H] and Grad(Age). 

The best linear fit to the Grad(Age) versus Grad[M/H] 
trend produces a slope of —1.14 ± 0.30. This result is 
comparable with findings of Sanchez-Blazquez et al. (2007; 
— 1.33 ± 0.37). However, other works on early-type galaxy 
samples give somewhat different results, like Mehlert et 
al. (2003; -0.95±0.11) and Rowle et al. (2010; -0.69±0.10. 
The concordance with Sanchez-Blazquez et al. result might 
be attributed to the fact that their sample, like ours, has an 
important fraction of disky, low-cro galaxies. Although dif- 
ficult to interpret in detail, this correlation indicates that 



large differences in star formation history between the inner 
and outer parts of galaxies yield large differences in chemical 
enrichment, which in retrospect is perhaps not so surprising. 



5 SUMMARY 

Using deep VLT-F0RS2 spectroscopy we have derived 
absorption-line index gradients for a sample of SO galax- 
ies in the Fornax Cluster spanning a wide range in galaxy 
mass. Combining these observations with the latest spectral 
synthesis models (Vazdekis et al. 2010), we have studied the 
properties of the stellar populations in these galaxies as a 
function of galactocentric distance, reaching the outer disk- 
dominated regions. 

Our main conclusions are: 

(i) The age gradients are positive or fiat, indicating that 
the star formation ended first in the disk-dominated outer 
regions of the galaxies. This centrally-concentrated star for- 
mation at later times could serve to enhances the bulge- 
to-disk ratios of these systems, helping the morphological 
transformation of spirals into SOs. 

(ii) The cessation of star formation in the central regions 
happened at earlier times in the brightest galaxies, consid- 
erably reducing their H/3 and age gradients. 

(iii) The age in the outer disk-dominated regions is > 
10 Gyr for the brightest and most massive galaxies and 
~ 4-5 Gyr for the least massive ones. Thus, star formation 
stopped first in the most massive galaxies not only in their 
centres but also in their oldest outermost regions, suggesting 
that "downsizing" is a rather universal phenomenon. 

(iv) Metallicity gradients, when detected, are always neg- 
ative, such that the galaxy centres are more metal-rich. This 
result fits with a picture in which star formation continued 
in the central regions, with enriched material, after it had 
stopped in the outskirts. 

(v) For the brightest galaxies, gradients in age and metal- 
licity are anti-correlated. This trend cannot be attributed to 
correlated errors. It indicates that large differences in star 
formation history between the inner and outer parts of SO 
galaxies yield large differences in chemical enrichment. 

In combination with the previous results obtained in 
this study (Papers I, II and III), a coherent picture is emerg- 
ing in which SO galaxies can be interpreted as spiral galaxies 
in which star formation has ceased. This conversion process 
seems to have occurred earliest in the most massive galax- 
ies as a further manifestation of downsizing. What we have 
been able to add by studying the radial gradients in these 
systems is the evidence that this cessation happened first 
in the disks of these systems, with star formation shutting 
down from the outside inward. Such a sequence fits well 
with the requirement for a process that enhances the cen- 
tral bulges of SOs relative to their disks when compared to 
unquenched spiral systems. 
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Figure 4. Gradients in Age, [M/H] and a-enhancement tracer versus each other and Central Age, [M/H] and a-enhancement tracer 
versus each other. Symbols as in Fig.[T] Continuous lines correspond to the best linear fits. Dashed lines correspond to flat gradients. In 
Grad{Age) vs. Grad[M/H] (top-left panel), errors as 1 cr contours have been drown for some galaxies. 
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APPENDIX A: LINE INDEX GRADIENTS 

In Fie:s. lAl| - |A5l radial measurements (along the semi-major 
axes) for 10 line indices are presented for the galaxy sam- 
ple. The radii are scaled to the total galaxy half-light radius 
(Rhalf) as calculated from Ks-band photometry. The radial 
profiles were folded around the galaxy centre (open and filled 
symbols; see each figure for details) . Linear fits are shown for 
the radial index profiles (indices in magnitudes) with those 
points lying in the central region affected by seeing excluded 
from the fit. The vertical lines show different radial scales 
for each galaxy when available. The dashed lines (in green) 
mark the bulge effective radius (IRe). The vertical contin- 
uous lines (in red) mark the disk scale- length (IRd). The 
dot-dashed lines (in blue) correspond to the radius (Rlim) 
from which the light profile is dominated by the disk com- 
ponent as calculated in Paper I. When Rlim is accompanied 
by horizontal arrows, it indicates that the line position is 
indicative only as its real location is too close to the galaxy 
centre to be shown. The resulting values of the index gradi- 
ents and their uncertainties are presented in Table IXT] 
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Figure Al. Ten index gradients along the semi-major axis of NGC 1316 and NGC 1380. Filled (open) symbols correspond to; SW (NE) 
(for NGC 1316) and S (N) (for NGC 1380). Central bin indices are not shown. Linear fits (error- weighted gradients) are shown for the 
radial index profiles where points within seeing-size from the nucleus where not used. Points with "*" symbols in Mg2 and Mgb indices 
where not used in the fit as they lie in a bad-pixel area of the detector. See the Appendix text for a description of the vertical lines. 
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Figure A2. Same as Fig. |XT]for NGC1381 and NGC 1380A. Filled (open) symbols correspond to; SE (NW) (for NGC 1381) and S 
(N) (for NGC 1380A). Points with "*" symbols in Mg2 and Mgb indices where not used in the fit as they lie in a bad- pixel area of the 
detector. 
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Figure A3. Same as Fig. lAll for NGC 1375 and IC 1963. Filled (open) symbols correspond to: tlie W (E) profile projected in the sky 
(for both galaxies). 
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Figure A4. Same as Fig.[XT]for ESO 358-G059 and ESO 358-G006. Filled (open) symbols correspond to the SE (NW) (for ESO 358-G059) 
and SW (NE) (for ESO 358-G006). 




Figure A5. Same as Fig. lAll for ESO359-G002. Filled (open) symbols correspond to the SW (NE) profile projected in the sky. 



Table Al. Line Index Gradients of SO galaxies in Fornax Cluster. 



Name 


Grad(H/3*) 


Grad(Fe5015*) 


Grad(Mgi) 


Grad(Mg2) 


Grad(Mg6*) 


Grad(Fe5270*) 


Grad(Fe5335*) 


Grad(Fe5406*) Grad(Fe5709* ) Grad(Fe5782*) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


NGC 1316 


-0.012 (0.007) 


0.030 (0.010) 


-0.008 (0.002) 


0.036 (0.003) 


-0.048 (0.008) 


-0.049 (0.011) 


-0.037 (0.016) 


-0.011 (0.011) 


-0.004 (0.004) 


-0.002 (0.010) 


NGC 1380 


0.009 (0.005) 


-0.011 (0.008) 


-0.026 (0.002) 


-0.043 (0.003) 


-0.060 (0.007) 


-0.013 (0.008) 


0.001 (0.013) 


-0.012 (0.009) 


-0.004 (0.003) 


-0.018 (0.009) 


NGC 1381 


0.005 (0.005) 


0.008 (0.005) 


-0.016 (0.006) 


-0.025 (0.002) 


-0.045 (0.005) 


-0.012 (0.007) 


0.003 (0.011) 


0.010 (0.008) 


-0.002 (0.003) 


0.010 (0.007) 


NGC 1380A 


-0.026 (0.007) 


-0.014 (0.009) 


-0.004 (0.002) 


-0.006 (0.003) 


-0.004 (0.007) 


-0.011 (0.010) 


-0.008 (0.015) 


-0.010 (0.011) 


0.003 (0.004) 


-0.008 (0.010) 


NGC 1375 


-0.061 (0.010) 


-0.058 (0.012) 


-0.001 (0.003) 


-0.020 (0.004) 


0.005 (0.010) 


-0.042 (0.014) 


-0.075 (0.022) 


-0.046 (0.016) 


0.005 (0.006) 


-0.019 (0.015) 


IC 1963 


—0.013 (0.006) 


0.015 (0.009) 


0.000 (0.005) 


—0.002 (0.003) 


—0.001 (0.007) 


0.006 (0.010) 


0.008 (0.016) 


0.011 (0.010) 


0.002 (0.004) 


-0.001 (0.009) 


ESO368-G059 


-0.020 (0.035) 


0.013 (0.047) 


0.018 (0.014) 


-0.024 (0.013) 


-0.073 (0.037) 


-0.026 (0.051) 


-0.086 (0.081) 


0.020 (0.058) 


0.009 (0.021) 


-0.047 (0.050) 


ESO368-G006 


0.003 (0.011) 


0.014 (0.019) 


0.047 (0.000) 


0.018 (0.070) 


0.000 (0.013) 


0.012 (0.018) 


-0.011 (0.028) 


0.006 (0.019) 


0.005 (0.007) 


-0.003 (0.017) 


ESO369-G002 


-0.138 (0.031) 


-0.021 (0.039) 


-0.015 (0.010) 


0.044 (0.012) 


0.097 (0.032) 


0.104 (0.046) 


-0.018 (0.073) 


0.005 (0.048) 


0.006 (0.019) 


-0.029 (0.049) 



Co 

a 

Notes: From (2) to (11), 1 cr RMS errors between () ; Col (1), galaxy name; Col (2), H/3 index; Col (3), Fe5015 index; Col (4), Mgi index in magnitudes; Col (5), Mg2 index in o 
magnitudes; Col (6), Mg6 index; Col (7), Fe5270 index; Col (8), Fc5335 index; Col (9), Fe5406 index; Col (10), Fe5709 index; Col (11), Fe5782 index. Gradients ealeulated with line "| 
indices in magnitudes (Paper III) with respect to log(R/RHALF) (radii normalized by the half-light radius of each galaxy). Significant gradients (^ 3(T detections) are highlighted in 
bold. 

o 

a 

■Si . 

;3 

■Si . 

a 

Co 
C2 

Co 

h— ' 
Or 



16 A.G. Bedregal, N. Cardiel, A. Arag on- Salamanca and M.R. Merrifield 



APPENDIX B: SSP PARAMETERS 

Figs.[Bl] - |B5] show SSP ages and [M/H] derived us- 
ing Vazdekis et al. (2010) models, and the a-element 
abundance tracer (Mg&/{Fe) ratio). The notation 
"log(Age[yr]) ([MgFe]' - H/3)" and "[M/H] ([MgFe]' 
- H/3)" is used to make explicit the fact that age and [M/H] 
have been derived using [MgFe]' - H/3 diagrams such as the 
one in Fig. (2] left panel. The solid lines (in black) show the 
best linear fits to the radial profiles. The slopes of these fits 
and corresponding errors are given in Table iBll 

The horizontal dashed lines in [M/H] plots represent so- 
lar metallicity. The two horizontal dashed lines in the [a/Fe]- 
tracer plots delimit the log(Mg6/(Fe)) range corresponding 
to [Q:/Fe]=0 [for metallicities —1.35 ^ [Z/H] ^ 0.35 and ages 
between 3 Gyr ^ t 15 Gyr according to Thomas, Maraston 
& Bender (2003) models; for NGC 1380 and 1381 the range 
is limited to models of ages between 8 Gyr ^ t ^ 15 Gyr]. 
Horizontal lines as in Figs. lAll - IA5l 

This paper has been typeset from a T^j]X/ ffl^jX file prepared 
by the author. 
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Figure Bl. Gradients along tlie semi-major axis of NGC 1316 (top) and NGC 1380 (bottom), including predicted Ages, Metallicities 
and a-element tracer values. Filled (open) symbols correspond to: SW (NE) (for NGC 1316) and S (N) (for NGC 1380). The description 
of the different lines can be found in the text. 
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Figure B3. Same as Fig. lBll for NGC 1375 (top) and IC 1963 (bottom). Filled (open) symbols correspond to: the W (E) profile projected 
in the sky (for both galaxies). 
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Figure B4. Same as Fig. [Bl] for ESO 358-G059 (top) and ESO358-G006 (bottom). Filled (open) symbols correspond to the SE (NW) 
(for ESO358-G059) and SW (NE) (for ESO 358-G006). 
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Figure B5. Same as Fig. IBll for ESO359-G002. Filled (open) symbols correspond to the SW (NE) profile projected in the sky. 
Table Bl. Age, [M/H] and Mgfe/(Fe> Gradients in Fornax SOs 



Name 
(1) 


Grad(logAge[yr]) 
(2) 


Grad[M/H] 
(3) 


Grad(«) 
(4) 


NGC1316 




"•""-0.10 


-0.02 (0.02) 


NGC 1380 


-0.08in! 


-0 21+"'^** 
"■^^-0.03 


-0.08 (0.02) 


NGC 1381 


-0 03+°"'' 


'^■^°-0.04 


-0.07 (0.01) 


NGC 1380A 


'^■^'^-0.08 


-0 34+0 14 
'^■•^^-0.11 


0.01 (0.02) 


NGC 1375 


46+'''3^ 
'^■^°-0.04 


n 72+0.12 


0.10 (0.03) 


IC 1963 


"■'-'--0.12 


-0.261°-^ 


-0.01 (0.02) 


ESQ 358-G059 


'^•■^'^-0.64 


-0.68tl;« 


-0.11 (0.11) 


ESQ 358-G006 


'-'■'^'-0.20 


0.08t«jl 


-0.00 (0.04) 


ESQ 359-G002 


QQ2 + 0.41 

'-'•y^_0.54 


_Q4. + 1.23 


0.16 (0.13) 



Notes: For (2) and (3), asymmetric la RMS uncertainties included; For (4), la RMS errors between parentheses; Col (2) and (3): SSP 
Age and [M/H] from H/3 vs. [MgFe]' diagram; Col (4): Gradient in log(Mgb/(Fe)) as an a-abundance tracer. Significant gradients 3 a 
detections) are highlighted in bold. 



